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Abstract 



A variational framework is developed to examine the equilibrium states of a semiflexible polymer 
that is constrained to lie on a fixed surface. As an application the confinement of a closed polymer 
loop of fixed length 2tvR within a spherical cavity of smaller radius, i?n, is considered. It is shown 
that an infinite number of distinct periodic completely attached equilibrium states exist, labeled 
by two integers: n — 2,3,4,... and p = 1,2,3, the number of periods of the polar and 
azimuthal angles respectively. Small loops oscillate about a geodesic circle: n = 2, p = 1 is the 
stable ground state; states with higher n exhibit instabilities. If R > 2i?o new states appear as 
oscillations about a doubly covered geodesic circle; the state n = 3,p = 2 replaces the two-fold as 
the ground state in a finite band of values of R. With increasing R, loop states make a transition 
from oscillatory and orbital behavior on crossing the poles, returning to oscillation upon collapse 
to a multiple cover of a geodesic circle (signaled, respectively, by an increase in p and an increase 
in n). The force transmitted to the surface does not increase monotonically with loop size, but 
does asymptotically. It behaves discontinuously where n changes. The contribution to energy 
from geodesic curvature is bounded. In large loops, the energy becomes dominated by a state 
independent contribution proportional to the loop size; the energy gap between the ground state 
and excited states disappears. 

1 Introduction 

Understanding how surfaces may constrain the configuration of biopolymers on mesoscopic scales is 
important in a number of processes in cell biology. It is particularly relevant in the packing of DNA 
within viral envelopes [U |2]. Modeling all of the relevant interactions is complicated: one needs to 
accommodate the competition between polymer elasticity and entropy; in the case of DNA, electric 
fields will be important. A nice short review in this context is provided in |3J. Recently, simulations 
treating various features of the confinement process have been performed [4, 5j; the former focuses on 
entropy, the latter on the dramatic effects of friction and the finite transverse dimensions of the confined 
object. 

In this paper, we address an aspect of the problem of a fundamental nature that does not appear to 
have been addressed previously in any detail: how does one characterize the equilibrium states of the 
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three-dimensional elastic energy of a polymer confined within a fixed surface? While this description of 
confinement leaves out a lot of the physics that is relevant in biological systems, it presents a well-defined 
problem exhibiting a striking level of complexity that is worth studying in its own right. 

The semi-flexible polymer will be modeled as a curve in three-dimensional space parameterized by arc- 
length s. The bending energy associated with a given conformation is quadratic in the Frenet-Serret 
curvature along the curve, «(s), 

H= 1 -JdsK(sf. (1) 

Curves of fixed length minimizing the unconstrained energy ([I]) were first studied in depth by Euler. 
A historical review is provided in Ref. [6J; a more contemporary approach to the problem is presented 
in [71 [8] and reviewed in Ref. [9J. An alternative framework, which lends itself to adaptation to the 
confinement problem, was developed in Ref. [10J. 

In the absence of constraints on a closed loop, there is not a lot to say: the only stable equilibrium is 
a circle. Confinement within a surface, however, will generally oblige the loop to adopt a non-circular 
shape, increasing its elastic energy. The contact with the surface itself may be complete or it may be 
partial; and contrary to one's initial guess, the bound state will not generally follow a surface geodesic; 
nor need it be unique. 

The wrapping of a semi-flexible polymer around a cylinder, a closely related problem relevant in the 
the winding of DNA around histone octamers, was first examined some time ago by Nickerson and 
Manning [TH EjQ Their strategy was to look at the independent degrees of freedom of the surface- 
bound polymer. While focusing directly on these degrees of freedom makes sense, it does not exploit 
the symmetries of the problem. For, even though the constraint breaks the Euclidean invariance of the 
three-dimensional bending energy, how this occurs is not arbitrary. A variational framework is developed 
here that involves the unconstrained degrees of freedom, imposing the constraint using a local Lagrange 
multiplier. This multiplier will quantify the loss of Euclidean invariance in the constrained system. Its 
value at any point along the loop will be identified as the local normal force that is being transmitted 
to the surface. 

The well-known integrability of the Euler-Lagrange equations for the unconstrained curve is a conse- 
quence of the Euclidean invariance of its energy. The constrained counterpart generally will not be 
integrable. In various interesting cases, however, the confining geometry will respect some subgroup of 
the Euclidean group. In particular, the conservation of torque associated with the residual rotational 
invariance of a sphere permits the Euler-Lagrange equation to be cast as a quadrature in the geodesic 
curvature which can be integrated directly. It also provides a direct recipe for the reconstruction of the 
loop from its curvature data adapted to the conserved torque vector as an axis of symmetry. 

In particular, we apply this framework to examine the confinement of a closed polymer loop of fixed 
length 2ttR within a sphere of radius R G < R. 

In contrast to an open polymer which will wind around a geodesic circle on the sphere when its length 
exceeds 7rR 0l the closure of the loop is incompatible with a geodesic unless R is tuned to be commensu- 

1 See also [13] and the work of Rudnick and Zandi [14]. A review is provided in Ref. [15]. There has also been 
some nice work done more recently by Van der Heijden et al. [16]. More directly relevant is the study of confinement of 
cylindrical sheets within a circular cylinder by Boue et al. [17]. 
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Figure 1: (Color online) State with two-fold symmetry for values of R in the interval [1,3] within a unit sphere: 
(a) and (b) display increasingly large oscillations about the equator; (c) oscillations develop overhangs; (d) first 
self-contacts made at the two poles; (e) self-intersecting triple orbit of sphere; (f) orbit collapses to triple cover 
of the equator. The normalized local confining force A is color coded in these figures. 

rate with an integer multiple of R . We show that there exists an infinite number of distinct completely 
attached states, labeled by a pair of integers, n and p\ n — 2, 3, 4 . . . , p — 1, 2, 3, . . . , the number of 
periods of the polar angle and azimuthal angles in one circuit of the loop respectively, n characterizes 
the dihedral symmetry with respect to the axis of symmetry; p characterizes the number of revolutions 
about this axis. Small loops oscillate symmetrically about a geodesic circle with p = 1 and an n-fold 

symmetry, n = 2,3, The two-fold is the stable ground state. For any finite values of R, the higher 

n-folds are unstable with respect to decay toward the two-fold. States with n = 2 are illustrated in Figs. 
[T](a) -[l](c) for various normalized values of R. 

Beyond some critical size, the loop will exhibit self-intersections on the sphere. When n = 2 this will 
occur when R = 2.127i? where the loop crosses the poles as illustrated in Fig. [ijd). We suppose that 
self-intersections are consistent with the physics and do not cost energy. 

As R is increased to R = 3R 0l the loop will collapse onto a geodesic circle which it will cover three times, 
(see Figs. [11(e) and [11(f). The discontinuity in the number of revolutions occurs at an intermediate value 
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of R where the loop crosses the poles. This is accompanied by a transition from oscillatory to orbital 
behavior. As R is increased above 3R 0l the dihedral symmetry of the lowest energy state descending 
from the n = 2 ground state will jump to n = 4 with the reestablishment of oscillatory behavior. 

The bending energy of this state as well as the total force that gets transmitted to the surface will be 
determined as a function of R. The three- dimensional bending energy decomposes on the surface into 
a sum of two terms: one is associated with the geodesic curvature, intrinsic to the sphere; the other is 
associated with the normal curvature inherited from the surface, constant on a sphere. The latter energy 
thus counts the number of times the loop is wound within the sphere and it is proportional to R. The 
geodesic energy periodically falls to zero whenever R = pR and the loop is geodesic. Local maxima, 
associated with the incommensurability of the loop length with geodesic behavior, are displayed between 
these values; we show that their values decrease monotonically with loop size. For large values of R, the 
normal contribution is always dominant. It does not depend on the state in question. Thus, in this limit, 
the energy gap between the ground state and excited states disappears. The local force transmitted by 
the completely attached ground state loop is positive everywhere. Surprisingly, the total force does not 
grow monotonically with loop length except asymptotically, where it grows linearly with R and coincides 
with the naive expression, energy divided by R . The change of symmetry as R passes through odd 
integral multiples of R 0l manifests itself in a positive jump in the transmitted force analogous to a Euler 
instability associated with the buckling into oscillations. 

When R > 2R 0l a new set of states appears as oscillations about a doubly covered geodesic circle, p = 2 
with an n-fold symmetry, n = 3,4, 5, . . . ; among these new states, the lowest energy is displayed in the 
three-fold. Also, because the geodesic energy is small when the oscillations are small, these states will 
have lower energy than the two-fold ground state in a finite band of values of R beginning at 2R . There 
will be similar behavior in an infinite set of bands of values of R where the descendant states become 
geodesic. One thus needs to reassess the stability of the two-fold ground state and its descendants. 
We provide a heuristic argument for stability by constructing a homotopy that interpolates between the 
states with n — 2, p — 1 and n = 3, p = 2. By examining the energy along this homotopy, we show that 
a steep energy barrier separates the two equilibrium states. The n — 2, p = 1 state and its descendants 
thus appear to be stable classically. For any R > 2R , there are two stable states that alternate as the 
ground state as the length is increased. 

The paper is organized as follows: In section [2j we describe the framework. In particular, the breaking 
of translational invariance will be quantified by the non-conservation of a vector along the loop. It 
will be shown that the Euler-Lagrange equation for the curve can be cast as the vanishing of a linear 
combination of the unconstrained Euler-Lagrange derivatives. The constraining force is identified as 
some other linear combination of these derivatives. In Sect. [3] the confinement of a closed loop by 
spheres will be considered. In Sees. [4] and [5] we analyze loops in the harmonic and non-linear regimes 
respectively. The equilibrium states of the loop will be identified and the forces they transmit to the 
surface determined. An assessment of the stability of these states will be provided. We conclude with 
a discussion and a few suggestions for future work in Sec. [6} A number of useful definitions, identities 
and derivations are collected in a set of appendixes. 
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2 Curves constrained to surfaces 



Consider a space curve r : s —> Y(s) parametrized by arclength constrained to lie on a surface E. This 
surface is described in parametric form by the mapping £ : (u 1 ^ 2 ) —> X(^ 1 ,^ 2 ). The confined curve 
can then also be described as a surface curve r s : s —> (E/ 1 ^), U 2 (s)). In order to enforce the condition 
that T lie on £, one adds to the energy H[Y] given by Eq. ([I]) a term enforcing this constraint: 

H C [Y, U a ] = H[Y] + Jds X(s) ■ [Y(s) - X([7 a ( S ))] , (2) 

where A is a vector-valued Lagrange multiplier defined along the curve. This constraint will break the 
manifest translational invariance of the energy H. 

The variation of H c with respect to the embedding functions Y can be cast in the form 

5 Y H C = J ds (F' + A) • 5Y , (3) 

where prime represents derivation with respect to arclength and the tension in the loop F is given by 

^k 2 -cJt + k'N + ktB . (4) 

Here {T,N,B} is the standard Frenet-Serret frame carried by the curve and r is its torsion. The 
constant c is associated with the constraint of fixed length. A derivation of Eq. ([?]) is provided in 



Appendix A 



In equilibrium, one finds that F' = —A. Thus, in the presence of the constraint, the tension in the loop 
is not conserved; the multiplier is identified as the external force associated with the constraint [IB] . 

The corresponding variation of H c with respect to U a (s) is given by 

SuH c = - [ ds\.e a 5U a , (5) 



where e a , a = 1,2 are the two tangent vectors to the surface adapted to the parametrization by u a . 
In equilibrium, A • e a = ; in equilibrium, the force on the curve associated with the constraint always 
acts orthogonally to the surface. Let us write A = An, where n is the unit vector normal to E. The 
combined result is that 

F' = —An . (6) 



An integrability condition for closed curves follows from Eq. ([6]): 

ds A n = . (7) 
This identity holds whether or not contact is complete. 

Using Eq. ([4]) along with the Frenet-Serret equations, a straightforward calculation decomposes F' along 
the normals 

F' = £ N N + e B B , (8) 



5 



where the Euler-Lagrange derivatives of the bending energy £ N and e& are given by 

£n = k" + k (y - t 2 -cj , (9a) 
SB = - (k 2 t)' . (9b) 

AC 

The tangential Euler-Lagrange derivative e T = F' • T vanishes identically, a consequence of the fact 
that the only relevant degrees of freedom are geometrical, whether the curve is constrained or not. 

The surface-bound curve also carries a Darboux frame, {T,n, 1 = T x n}. Relevant properties of this 



frame are summarized in Appendix B The Frenet normals are related to their Darboux counterparts by 



a rotation about the tangent vector: 

N = cos uj n + sin uj 1 ; B = — sin uj n + cos uj 1 . (10) 

The relationship T' = kN permits one to express the angle of rotation in terms of a curvature ratio: 
since; = k 9 /k or cos uj = K, n /K, where k q and K n are, respectively, the geodesic and normal curvatures 
along T E , 

fy = r-l, « n = T , -n. (11) 
Using these expressions it is possible to decompose F' given by Eq. ([8]) in a form adapted to the surface, 

F' = (coscj £n — sinoj £&) n + (since; £n + cos uj Sb) 1 ■ (12) 

The projection of Eq. ([6]) onto 1 provides the Euler-Lagrange equation in the remarkably simple form, 

£\ = sin uj £n H~ cos ujs& = . (13) 

The multiplier A does not appear. The corresponding projection onto n determines A, 

A = — cos uj £n H~ sin uj e& . (14) 



In this approach, one sees explicitly how both the Euler-Lagrange equation (13) and the confining force 



(14) are constructed out of the two unconstrained Euler-Lagrange derivatives. The normal force is 



completely determined when the local geometry is known. 

The Euler-Lagrange equations for an unconstrained elastic curve, given by £ N = and s& = 0, are 
replaced by the single equation, s\ = 0. The apparent discrepancy in the number of equations reflects 
the fact that a space curve possesses two independent modes of deformation whereas a surface curve has 
only one. In general, the integrability of the former pair of equations is surrendered when the constraint 
is present. 

2.1 Euler-Lagrange equation in terms of surface curvatures 



Using the identities ( |B.3| ) and (B.4) the Euler-Lagrange equation (13) can be expressed completely in 
terms of surface curvatures, 
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This agrees with the equation derived in Ref. [12J using a very different approach. Note that it involves 
the curvatures as well as the geodesic torsion and, in general, the curve will not follow a geodesic with 
Kg = 0. For a geodesic to minimize bending energy, one requires one of the following to occur: \<? n r g 
is constant; the curve coincides with an asymptotic line with n n = or a principal curve with r g = 0. 
Such conditions typically do not occur unless they do so trivially. 

The magnitude of the force A transmitted to the surface, given by Eq. ( [I4| ), assumes the form 

-» = <K(^-f t )+^. (16) 

Its magnitude will vary along the contact region, even for confinement by a sphere. This expression is 
missing in the framework presented in [12J. A curious consequence of the symmetric decomposition of 
the Frenet-Serret curvature in terms of the geodesic and the normal curvatures is the fact that A turns 



out to be identical to the Euler-Lagrange derivative e\ given by (15) under the interchange of k q with 
K n and a change of sign of the last term. 

2.2 Confinement and the loss of rotational invariance 

In general, under confinement, one surrenders not only translational invariance but also rotational in- 
variance. The torque about the origin per unit length of the curve, M, is given by 

M= YxF + S, (17) 



where S = —kB = K g n — K n l . The first term on the right in Eq. (17) is the torque due to the force 
F; the second term is the bending moment originating in second derivatives in the bending energy. For 
a free curve, M 7 = 0, which can be cast in the manifestly translationally invariant form S 7 + T x F = 
[18] . In general, M is not conserved. One has instead 

M 7 = ^(Yxl)-A(Yxn) . (18) 

Thus, in a confined equilibrium with e\ — 0, M will not generally be conserved. The source is given by 
the moment of the force associated with the constraint. 

If the confining geometry is a sphere centered on the origin, so that Y is directed along the normal 
vector, M will be conserved. If it is symmetric about some axis (say the z axis), then the conserved 
quantity is the corresponding projection of M, that is, M 3 = M • z. 



3 Spherical confinement 

Consider a closed curve of length S = 2ttR, confined within a sphere of radius R G . We normalize lengths 
in terms of R . On a sphere, the extrinsic curvature tensor is proportional to the metric K ab = g abl so 
that the normal curvature is constant, n n = —1, and the geodesic torsion vanishes, r g = 0. Using Eq. 
(A. 6), the vector F given by Eq. [4| then reduces to 



\_ . 1 



F= [-£+<T]T + K' g l where a = --c. (19) 
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F is everywhere tangent to the surface. It is also completely determined by the intrinsic geometry. 
However, it is not conserved. 



The torque vector defined by Eq. (fl7|) is given by 



M 



^ + (7-l)l + « fl n. 



(20) 



The rotational invariance of bending energy confined to a sphere implies that M is a constant vector. In 
particular, its length is a constant. This provides a quadrature for the geodesic curvature (M 2 = M-M): 



ft: 



M 2 = (k') 2 + + + k 



(21) 



It is simple to check that the condition (M 2 )' = in Eq. (21) reproduces the Euler-Lagrange Eq. (15) 
on a sphere: 



Kg + Kg 



KZ 



+ a 



0. 



(22) 



This equation is also identical to the one describing an elastic loop on a sphere with an energy density, 
ft 2 /2. The minimization of the constrained three-dimensional bending energy coincides in this case with 
that of two-dimensional spherical bending energy. This is a consequence of the decomposition of the 
Frenet-Serret curvature into normal and geodesic parts, k 2 = K? g + K 2 n as well as the fact that K n is a 
constant. Spheres are special in this respect. 

Various mathematical properties of elastic curves on spheres were described by Langer and Singer in the 
1980s [7]; see also [19] for a numerical treatment of the problem. Recently they have been reexamined 
in some detail in [20]. The connection to the description of conical defects in unstretchable flat sheets 
was developed in |21j and |22j. While the mathematical literature provides a useful point of departure, 
the absence of any reference either to energy or to the transmitted forces limits its usefulness. 



Note that in a geodesic, k 9 = everywhere. Equation (21) then implies that M 2 = (a — l) 2 . Unless 



R is an integer, however, geodesies will be inconsistent with the boundary conditions associated with 
closure. If k 9 = is accessible anywhere along the loop, there will be a non-trivial bound on M from 
below: 



(a - l) 2 < M 2 
In equilibrium, all confined loops will satisfy this bound. 



(23) 



It is also useful to cast Eq. (21) in the alternative form 



w 



here 



V(Kg) = 

and E 2 = M 2 + 2a — 1 is manifestly positive^ 



K„ 



(24) 
(25) 



-This is a weaker bound than Eq. (23) 
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If Kg is identified as the position of a particle and s as time, then Eq. (24) described the motion of 
a particle of mass m = 2 and total "energy" E 2 in the symmetric quartic potential V(k 9 ). If a > 0, 
V possesses a single minimum at k q — 0; if a < it possesses two symmetric wells centered at 
^ = 2a, separated by a local maximum at k q = 0. 

While the particle analogy is useful, it does have its limitations. E 2 also is not the energy of the loop 
and the potential depends on the constant of integration a which we are not free to tune but, like the 
'energy', is itself determined by the boundary conditions. 

The qualitative behavior of the loop will depend on the turning points of the potential, and thus on the 
relative values of a 2 and E 2 \ 



(1) E 2 > a 2 (equivalently Eq. (23) is satisfied). In this parameter regime, there are only two turning 
points. One has 

E 2 - V{k 9 ) = \(n* g + K 2 ) (k 2 - kJ) , (26) 

where K 2 = 2(E + a) > and k\ = 2(E - a) > 0. Thus rf g = when k q = ±k x . The geodesic 
curvature thus ranges in the symmetric interval [— fci,fci]. This will be independent of the sign of a. The 
loop will oscillate symmetrically about the equator where k 9 = (see Eq. ([29])) so that § ds^g = 0. 

(2) E 2 < a 2 . This regime is inaccessible physically if the loop is closedj^] 

The reconstruction of the loop from its curvature data involves examining the conserved torque vector 
M. Without loss of generality, it is always possible to allign M along the z axis, M = Mz (we follow 
Ref. |22j). The normal vector n is parametrized in terms of spherical polar coordinates $ and tp t 

n(s) = (sin i?(s) cos (p(s), sin i?(s) sin ip(s), cos i?(s)) . (28) 

The projection of M, given by Eq. (20), onto n determines the polar angle in terms of 

M-n = Mcostf = /^; (29) 

its projection onto 1 determines the azimuthal angle 

M-l = -Msin 2 V = - (— • (30) 



Thus the projections of M onto the Darboux frame determine the embedding functions of the curve on 
the sphere in terms of k q and the two constants, a and M. 

ip will not generally increase monotonically along the loop. It will exhibit overhangs with tp f = 0, if 
a < 1. Combining Eqs. (29) and (30) gives 

M / M 2 + 2(q-l) \ 
^ = TT ~2 1 • ( 31 ) 



2 V M 2 -K 2 g 



3 It can only arise if a < so that the potential possesses two wells and the trajectory in n g is confined to oscillate in 
one of them. One can write 

£ 2 - V( Kg ) = \ (4 - k]) (kI - nl) , (27) 

where k\ = — 2(E + o) > and n 2 _ = 2(E — a) > 0. k q is then confined to the lie in one of two intervals with a definite 
sign. It is thus confined to inhabit a single hemisphere. It is intuitively clear that any such state will spontaneously unbind 
into the interior of the sphere where it may relax into a lower energy bound state. 
4 The projection onto T = n' reproduces (the derivative of) Eq. (29). 
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Note that Eq. (29) implies the bound on K g , \k 9 \ < M, already implicit in the quadrature, Eq. (21)J^j 
As we will see, M will always be bounded. Thus k 9 and with it the geodesic energy will also be bounded. 



Equation (29) also implies that the extremal values of k 9 occur where the polar angle is turning. The 



bound is saturated when the loop passes through the poles. 



Using the identification, (29), it is possible to recast the quadrature in terms of 9. One has 



e ,2 + 



1 /M 2 + 2(a-l) 



M 



sin# 



M 2 sm6) = 1. 



(32) 



In this form, it is clear that access to the poles is possible only when M and a are tuned so that the 



coefficient of the divergent term in the potential appearing in Eq. (32) vanishes 



M 2 + 2(<j — 1) = . 



(33) 



This will require a < 1. 



When Eq. (33) is satisfied, the evolution of simplifies. Equation (31) assumes the form <p' = —M/2 
so that (f increases linearly with arclength along the loop. In general, the dependence of tp on s is not 
monotonic, much less linear. This identity, in turn, implies the value M = 2/R (and as a consequence 
of Eq. (33), a = 1 — 2/R 2 at pole crossing. 



Let us first suppose that the loop is sufficiently small so that neither self-contact nor self-intersections 
occur. The periodic motion in the potential V(k 9 ) implies an n-fold dihedral symmetry: the loop closes 
upon completing n periods of 9 in one revolution of the polar axis so that 9(s + 2nR/n) = 9{s) and 
(f(s + 27rR/n) = (p(s) + 2n /n, where n > 2 is an integer]^] Closure provides a quantization of physical 
states H 

In equilibrium, the number of circuits of the polar axis, p, will increase with loop length, so that the 
boundary condition on is replaced with Lp(s + 2'nR/n) = (p(s) + 2 r K r pjn. Using Eqs. (29) and (31), 



and the quadrature (24) it is possible to cast the boundary conditions in the form 

2nR f kl dK n 



An 



and 



2np M 



An 



ki 



V^ 2 - Wg) ' 

kI + 2(a - 1) 



2 J _ y£ 



M 2 



K„ 



(34) 



(35) 



where k\ is the turning point of the potential, defined by Eq. (26). Equation (34) is independent of 



p, whereas Eq. (35) is independent of R. Together, they determine the two constants of integration a 
and M in terms of the loop radius R and the two integers n and p. 



5 A sharper, if less transparent, bound follows from Eq. (24) which implies, for positive a, k 2 < 2^2+2(7-1-2(7. 
6 The identity (29) implies that the former is equivalent to n g (s + 2i\Rjn) — k 9 (s). The quadrature then implies that 
closure is smooth. 

7 A one-fold n = 1 is incompatible with the four-vertex theorem for a sphere so does not occur. 
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4 Weak confinement by Spheres 



While Eq. (21) can be integrated exactly in terms of elliptic functions [ZlElGQ]), a perturbative approach 
to the problem is instructive. Let us thus suppose that AR := R — 1 <C 1. Such a loop is sufficiently 
small that neither self-contact nor self-intersections occur. We thus expand the function k 9 as well as 
the constants a and M in powers of e = V AR, the small dimensionless parameter in the problem: 



Kg = tti + k 3 + • • • ; M = M + M 2 + • • • ; a = a + a 2 + • • • . (36) 

The equilibrium states are described by small oscillations about a geodesic circle on the sphere with 
Kg = 0. The harmonic approximation of the quadrature (21) about k q = reads 



(k[Y + a K( = M£ - ((7 - 1)' + 2M M 2 - 2a a 2 . 



(37) 



At lowest order, the arclength coincides with the azimuthal angle, p, and the geodesic curvature along 
a closed loop is given by 

K\{s) ~ A\ cos np , (38) 

where n is an integer and A\ is a constant. Ai? is determined by the amplitude A\ and the n (see 
Appendix C). 

The quadrature implies that a = n 2 and M = n 2 — 1. It also implies the constraint 

^(n 2 -l)(M 2 -a 2 ) = ^ (39) 
on the difference of their second order corrections. 



Equation (38) implements the boundary conditions at lowest order. To complete the specification of cr 2 



and M 2 in terms of AR, it is necessary to examine the boundary conditions (34) and (35) correct at 



next to leading order. One finds that, for n > 2, Eqs. (34) and (35) together imply 



and 



AR 



AR 



4n 4 



(<t 2 - M 2 ) 



4^ + A '«-l<^ 2 



A 2 



-(<x 2 -M 2 ) 



1 



n 2 + l i2 



4(n 2 - 1) n 2 - 1 



(40) 



(41) 



The details of the derivation are provided in Appendix C. Using Eq. (39) in (41) one reproduces the 
relationship Eq. (Oh between A l and AR, A\ = 4(n 2 - 1)AR. It then follows that, for n ^ 1, 



<x 2 = ^(3-7n 2 )Ai?, M 2 = ^(l-n 2 )Ai?. 



(42) 



4.1 Energy and transmitted force 

The bending energy of the loop confined by the sphere decomposes into a sum of geodesic and normal 
parts, reflecting the decomposition of the Frenet curvature, k 2 = k 2 + 1, 

H= l -jds(K 2 g + l) := H g + H n . (43) 
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H n = 7iR is the energy associated with an elastic rod that been wound into a circular coil of radius 
Ro = 1. It grows linearly with loop length 2ttR and is state independent. For a weakly confined loop 

H/H\ oop w l + 2n 2 Ai?, (44) 

where H\ oop = tt/R is the bending energy of a circular loop of radius R. The energy increases linearly 
with loop size. This will not be true in longer loops. 

For a fixed value of AR, the energy increases quadratically with n. The ground state, as one would 
have predicted, is the completely attached n = 2 state. As we show, all states with n > 3 are unstable. 



The force transmitted to the sphere at any point, given in Eq. (16), takes the particularly simple form 

A = i^ + a. (45) 



It differs from the local energy density: compare Eqs. (43) and (45). Its spatial dependence is completely 
determined by the geodesic curvature. It is bounded from below by a. This lends a direct physical 
interpretation of a. For small AR and each n > 2, the force is given by 

AR 

\ n « n 2 + — [3 - 7n 2 + 4 (n 2 - l) cos 2 nip] ; (46) 

while it oscillates along the loop, it remains positive everywhere. \ n does not vanish in the limit AR — >> 0. 
This is interpreted as an Euler instability. The buckling of the circular loop into an n-fold involves a 
critical compression in the loop provided by this normal force. Thereafter, A n decreases linearly with 
AR. This may appear counterintuitive. This behavior will be put in context when we examine loops 
outside of perturbation theory. 

4.2 Stability 

Here we examine the stability of the weakly confined n-folds that have been described. To do this 
we require the second variation of the energy with respect to small deformations of the loop. This is 
expressed in the form [23] 

5 2 H = [ dip$C<S>, (47) 



where $ is the deformation along 1, and the self-adjoint operator C is given by 



&A +(n 2 + l)^-+n 2 . (48) 



<V v J dip 1 

The fixed length constraint along the spherical surface implies a global constraint on the normal defor- 
mation, 

dsK g $ = . (49) 

Periodicity implies that the modes of deformation are represented by a constant (m = 0), smrwp 
and cosm^, m = 1,2,3, ... . For fixed n, the eigenvalues of C are then labeled by the integer m = 
0, 1,2, ... , given by 

C m = (m 2 -n 2 )(m 2 -1). (50) 
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The constant mode with m = has positive C = n 2 . All eigenvalues with m > 1 possess a two-fold 
degeneracy. 

There are four zero modes satisfying = 0; two occur at m = 1 and two at m = n. The mode 
siring oc corresponds to rotation of the loop about the axis of symmetry. The two modes with m = 1 
correspond to rotations about an orthogonal axis. These three modes are the zero modes anticipated by 
the rotational invariance of the bending energy. The fourth zero mode cosriLp is inconsistent with the 



fixed length constraint (49) and, so, is unphysical. It is also the only mode of deformation inconsistent 



with this constraint in this regime. 

The two-fold ground state with n = 2 is stable. There are no modes of deformation with negative 
eigenvalue. 

All excited confined states are unstable. There will be 2{n — 2) unstable modes of deformation cor- 
responding to m — 2, . . . , n — 1 lying between the zero modes at m — 1 and m = n. All modes of 
deformations with m > n contribute a positive energy. The first excited state with n = 3 has two modes 
of decay, of equal energy, into the ground state. The dominant mode of instability in higher energy 
states is not directly towards the ground state involving a cascade of instabilities. 



5 Strong spherical confinement 

Let us now examine the shape adopted by a loop of finite R confined by the sphere. 



Equation (22) can be integrated in terms of elliptic functions to give k 9 as a function of s [7, 9, 

K g (s) = ft cn [qs,m] , k q = 2\fmq. (51) 

The function cn[x,m] is the Jacobi elliptic cosine [24]. The angular wavenumber q is given in terms of 
the constant E defined below Eq. (25) by q = \[E\ the modulus m is defined by 



m = - [ 1 
2 



~n2 



(52) 



The curvature depends on the two parameters a and M through the parameters q and These 
parameters will be determined explicitly in terms of R and n using the boundary conditions associated 
with the closure of the loop. 

Using the fact that the period of cn is given by 4/C[m], where JC[m] is the complete elliptic integral of 
the first kind |24j, it is possible to cast the boundary condition on K gt given by Eq. (34), in the form 

q = 4n)C[m\/(27vR). (54) 



Integration of Eq. (30) gives 



/ n. M ( q 2 + 1 



2 V^(^ 2 - 1) 



n 



Amq 2 ~| \ 

" ( g 2 _ !)2 ' am l^ m J > m - s ) 



3 Definitions of q and M are inverted to give 

a = q 2 (1 - 2m) , 



M 



(q 2 - l) 2 +4ra<? 2 . 



(55) 



(53) 



The modulus m will lie in the interval [0, 1]; this bounds a in terms of q\ \<j\ < q 2 , changing sign when m = 1/2. 
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where U[rj, am[x, to], m] is the incomplete elliptic integral of the third kind and am[x,m] is the Jacobi 
amplitude [24j. Closure of the loop after p circuits of the polar axis, Lp(2nR) = 2np, then reads 



2^M R ( q 2 + 1 

/C[m] (g 2 - 1) 



4mg 

m 



2np . (56) 



Modulo Eq. (54), Eq. (56) determines m implicitly as a function of R. This completes the formal 



construction of the confined loops. 

5.1 Ground state n = 2 f p = 1 and its descendants 

The two-fold state with n = 2 is illustrated in Fig. [I] for values of R in the interval [1,3]. When 
AR is small, it consists of a small oscillation about a geodesic circle (see Fig. 1(a)) consistent with 
the perturbative description. As R is increased the oscillations about this circle initially increase in 
amplitude as more of the surface of the sphere is explored (Figs. 1(b) and 1(c)). In Fig. 1 (c) overhangs 
with ip' = occur on the increasingly crowded loop. At a critical value R 2 = 2.127, the loop makes 
self-contacts at the two poles (Fig. l(d))j^] 

For values of R above R 2l one needs to look more carefully at the boundary conditions. If, as we assume, 
the microscopic physics accommodates self-intersections on the surface and they occur without costing 
energy, the mathematical curve described by the elliptic function continues to represent the physical 
Ioop0 

There is a qualitative change in the behavior as R is increased through R 2 . On crossing the poles, the 
loop makes two additional revolutions about the polar axis so that p = 1 must be replaced by p = 3 in 



the boundary condition Eq. fl35| ). To see this, note that the loop will cross the poles four times, twice at 
each pole. A discontinuity of — 2tt is introduced in at each crossing. These discontinuities contribute 
to its period; thus, the period of 2tt describing a single revolution gets replaced by 2tt — 4(27r) = — Qtt. 
The orientation is also reversed J 17 ] As Fig. [I] illustrates, this discontinuity marks a transition from 
oscillatory to orbital behavior. 

As R is increased further the two regions bounded by the self-intersections grow (Fig. 1(e)). When 
R — 3 the loop degenerates into a triple cover of a geodesic circle (Fig. 1(f)). 

5.1.1 States as trajectories in parameter space 



Each equilibrium state can be represented as a point on the parameter space (a, M). As R increases, 
the state will follow a trajectory R —> (cr(R) : M(R)) labeled by the two integers n and p^\ 

9 For each n there will be a critical value R n where this occurs, whose magnitude increases with n. As pointed out in 
[22] , this value saturates. Beyond some critical loop length, all equilibrium states are self-intersecting. 

10 lf self-intersection is prohibited, the physics will be very different. It has been described in the conical context in Ref. 
|22j . and explored numerically in detail in Ref. |25j . As the loop crowds itself on the sphere, there will be a steep rise 
in the energy associated with this self-confinement. Beyond some point, partially attached loop conformations would be 
expected to become energetically favored. 

11 For an n-fold, the corresponding number of revolutions is p = 2n — 1. 

12 The boundary condition on (p, Eq. (35), does not involve R explicitly. It thus identifies the trajectories on the (<r, M) 
parameter space without locating the position along these trajectories. The latter is provided by Eq. (34). 
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To facilitate the interpretation of trajectories in parameter space, it is useful to first locate relevant 
geometrical landmarks. Various significant parameter curves are represented in Fig. |2) 



(1) The parabola II given by M 2 = (a — l) 2 (saturating the inequality (23)) provides a boundary on 
parameter spacej^] Points on the boundary describe geodesic loops, with n g = and n' g = 0, covering 
the equator an integer number of times. Points below the parabola do not describe confined loops. 

(2) Points on the straight line L, given by M 2 = 2(1 — a), describe loops that pass through the poles. 
Points above this line represent loops that display oscillatory behavior; those below it display orbital 
behavior. 

(3) Overhangs occur to the left of the vertical line, a = 1. 

The solid curve, Li in Fig. [^represents the trajectory of the state with n = 2 and p = 1 as R is increased 
in the interval [1,3]. The initial state with R = 1-a single geodesic circle-is represented by the point 
Pi = (4,3) lying on the boundary parabola II. In this interval, both a and M decrease monotonically 
with R. The trajectory terminates in the point Qi = (4/9,5/9) on the boundary parabola, where the 
loop degenerates into a triple covering of a geodesic circlep] 

More generally, whenever R — p, where p is an integer, loops may form p-fold coverings of a geodesic 
circle. As R is raised above R = p, it is possible to examine the state perturbatively in a manner 
analogous to that of small loops. Suppose that the loop oscillates with an n-fold symmetry, so that 
Kg ~ cosnLp/p. Consistency with R > p places a lower bound on n: n > p. Thus one must have 
n = p + l ? p + 2, . . . p*l The state with p + 1-fold dihedral symmetry will minimize the energy among 
these states. It is thus identified as the ground state. In particular, when R passes through R = 3, the 
symmetry of the ground state changes from n = 2 to n = 4 (see Fig. [3]). 



It is straightforward to determine the values of a and M in the ground state with n = p + 1 as R 
passes through R = p\ a(p) = (p+ I) 2 /p 2 , and M(p) = (2p + l)/_p 2 , representing a point lying on the 
boundary parabola II. 

As R is increased the amplitude of oscillation increases. At some point the loop will (re)cross the 
poles where the number of revolutions changes by two. As R = p + 2 is approached, the loop morphs 
continuously into a p + 2 cover of a geodesic circle. In this region, the ground state can also be treated 
perturbatively. Now k 9 w cos((p+ l)/(p + 2)(p). The corresponding end point in parameter space as 
R approaches R = p + 2 is given by a = (p+l) 2 /(p + 2) 2 , M = (2p + 3)/(p + 2) 2 , which also lies on 
the boundary parabola II. 

The trajectory in the parameter plane representing the n = 2, p = 1 ground state loop and its minimum 
energy descendants is illustrated in Fig. [2j It consists of a number of disconnected curve segments, L p , 

13 For graphical purposes we plot M 2 vs a; hence "parabola". 

14 Prior to crossing L, the trajectory L\ is approximated accurately by the straight line, M 2 = 9 (1 + 6/25 (a — 4)), 
with a slope 54/25, predicted in perturbation theory. Perturbation theory thus provides an accurate description of the 
problem outside its expected rang e of validity. 

15 Note that the analog of Eq. (C.3) for p revolutions implies 



+ (57) 



4M 2 \p 2 
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cr-1 

(a) Ground States n = 2, p = 1 and descendants 



(7-1 

(b) Ground States n — 3, p = 2 and descendants 



Figure 2: Trajectories in parameter space (a — 1, M 2 ). (a) States with n = 2 are represented by the solid curve 
L\ for i? in the interval [1,3]. This curve originates at the point P\ = (3,9) when R = 1 (geodesic circle) 
and terminates at the point Q\ = (—5/9,25/81) when R% (a triply covered geodesic circle). Both lie on the 
boundary parabola II. The straight line L indicates loops passing through the poles. The descendant states with 
n = 4 [n = 6) which occur for R in the interval [3,5] ([5,7]) are represented by the dashed curve L3 (dotted 
curve L5) in (a), where they are also shown with zoom in the inset. The counterpart of (a) for the state n = 3, 
p = 2 and its descendants is represented in (b). States with n = 3 for R in the interval [2,4], as well as the 
descendants with n = 5 and n = 7 for R in the intervals [4,6] and [6,8] respectively, are represented by the 
solid, dashed, and dotted curves, £2,^4, Lq. 

p = 1, 3, 5, . . . with end points P p = ((p + l) 2 /p 2 , (2p + l)/p 2 ), and Q p = ((p + lf/{p + 2) 2 , (2p + 
3)/(p + 2) 2 ) lying on the boundary parabola II. The points Q p and P p+2 describe loops with R = p as 
it is approached from below and above respectively. 

There are discontinuities in a and M 2 at i? = p = 3, 5, 7, . . . and the loop passes through a p-fold 
covering of a geodesic circle. These discontinuities are identified analytically using perturbation theory. 
Their origin is the transition from orbital back to oscillatory behavior with a change from p — 1-fold to 
p + 1-fold symmetry. 

As p -> oc, P p , Q p -> (1, 0) so that a -> 1 and M -> 0. 
In Fig. [4] a and M 2 are plotted as functions of R. 

5.2 Ground state n = 3, p = 2 and its descendants 

The n-fold symmetry states with p = 1 and their descendants do not exhaust all possible self-intersecting 
states. If R = 2, a set of states consisting of a double covering of a geodesic circle (p = 2) comes 
into existence with symmetry n = 3, 5, 7, . . . . The sequence L 2 with n = 3 is illustrated in Fig. [5] for 
values of R in the interval [2,4]. As R is increased, the state with n = 3 has least energy. As before, 
one can show that it is the only stable member among these states. It also has a lower energy than its 
counterpart with n = 2, and p = 1 illustrated in Fig. [TJ 

The trajectory L 2 representing these states in parameter space as R is varied in the interval [2,4], is 
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(a) R = 3.05 



(b) R = 4.06 



(c) R = 4.99 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.9 L 



Figure 3: (Color online) States with four-fold symmetry, n = 4, p = 3 for R in the interval [3, 5]: (a) increasingly 
deformed triple covering of the equator; (b) pole crossing, four points at a time); (c) state making five orbits 
of sphere before it becomes a quintuple covering of equator. As in Fig. 1 colors represent the normalized local 
confining force. 




1 2 3 4 5 6 7 1 2 3 4 5 6 7 

R R 
(a) (b) 



Figure 4: a (a) and M (b) vs. R for n = 2,p = 1 and its descendants in the interval [1, 7], as described in the 
text represented by the trajectories Li,Ls and L5 in Fig. [2] Asymptotically, a —> 1 and M — » 0. 



17 





(a) R = 2.05 




(b) R = 2.5 




(c) R = 3.0 




(d) R = 3.08 



(e) R = 3.5 



(f) R = 3.99 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



0.9 1. 



Figure 5: (Color online) State with three-fold symmetry for R in the interval [2,4]: (a),(b),(c) increasingly 
deformed double covering of the equator; (d) pole crossing, three points at a time); (e) state making four orbits 
of sphere; (f) quadruple covering of equator. As in Fig. 1 colors represent the normalized local confining force. 

represented in Fig. |2][b) by a solid curve. When R = 4 a transition occurs to a state with a five-fold 
symmetry. The sequence of trajectories L 2 , L 4 , L 6 , ... is generated. Their endpoints converge to the 
same point = (1, 0) as those of the sequence L 1? L 3 , L 5 , . . . . 



5.3 Loop Energy 



The bending energy Eq. (43) can be recast explicitly as a function of R: 

16 n 2 



H 



7lR 



K[m] (S[m\ + (m- l)K[m]) + ttR , 



(58) 



where S[m] is the complete elliptic integral of the second kind |24j, and m is determined by solving Eq. 



(56) 



H is plotted as a function of R for the ground state with n = 2 and p = 1 and its descendants in Fig. 
|6| Expanding expression (58) to second order in AR = R — 1 we reproduce Eq.(44). 
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4 

R 



Figure 6: The total energy of the n = 2 ground state and its descendants as a function of R is represented 
by a solid curve. It is sum of geodesic and normal parts which are plotted separately as a dashed curve and a 
dotted line. The normal energy grows linearly and dominates in large loops. By contrast, the geodesic energy 
is bounded, vanishing at odd integers where its derivative is discontinuous. Its maxima between zeros decrease 
inversely with R. 



The normal energy, represented by the last term, grows linearly with R and is state independent. In 



contrast the energy associated with geodesic curvature, represented by the first two terms in Eq.(58), 
depends sensitively on the state. It is bounded and its maxima fall monotonically as the loop becomes 
large. It is simple to place a crude upper bound on this falloff: one has M « 2/R, so that the 
H g < 2/R. Thus in a large loop, the geodesic contribution to the energy is negligible compared to its 
state independent normal counterpart. 

The geodesic energy vanishes when R = 1,3,5,..., where the loop collapses to a multiple covering 
of a geodesic circle. The discontinuities in the derivative of H with respect to R at R = 3, 5, . . . , 
are directly associated with the change of symmetry at these values of R. In the intervals between 
consecutive values, H increases linearly from zero and rises to a maximum value before falling to zero, 
again linearly]^] Behavior is not symmetrical in these intervals. The slopes are different at the two ends 
and the maximum is not centered. While the maxima are not exactly periodic, they do migrate to the 
center of their respective intervals as R increases. The existence of these maxima can be understood as 
a consequence of the incommensurability of the loop with geodesic behavior. The strongly asymmetrical 
initial maximum is associated with the development of overhangs on the loop before it crosses the poles. 
As the loop becomes longer this incommensurability plays a diminishing role. 

This qualitative behavior is repeated for the odd ground states n = 3, p = 2 and its descendants, as well 
as the excited counterparts of these states. Strikingly, the energy of any closed equilibrium loop tends 
to a common value in the limit, independent of the state, that completely dominated by the normal 
energy. The finite gap between the ground state and excited states state disappears. 

If R — 1 is not small, the stability analysis is rather more complicated than the one presented in the 
context of weak confinement. The results of Ref. [23] in another context implies that, although the 



16 The initial linear behavior was described in perturbation theory. The existence of a local maximum was not. 
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details differ, the conclusions for weak confinement continue to hold for values of R before the onset of 
self-intersection and only the state n — 2, 3, 4, . . . , p = 1 are excited. A detailed analysis of stability 
has yet to be performed beyond this point. One would, however, expect that the instabilities persist. 

The stability of the n = 2, p = 1 ground state and its descendants needs to be reassessed when R>2. 
For now one has to accommodate the existence of a set of states described by n = 3, 4, 5, . . . p = 2 and 
their descendants. In particular, within a finite interval of values of R above R = 2, the state n = 3, 
p = 2, represented in Fig. [5] has lower energy than the state with n = 2 and p = 1 represented in 
Figure [T]). See Fig. [7| The difference in energy is associated with the significant geodesic curvature of 
the latter when R « 2. Its descendants with n = 5, 7, . . . also will have lower energy than their even 
counterparts when R « 4, 6, 

There are, of course, no smooth deformations taking one from a loop with p — 1 to one with p = 2 
that remain attached to the sphere. A hairpin costing an infinite bending energy will necessarily always 
form. However, it is possible to sidestep this topological obstruction by permitting the loop to detach 
into the interior. While a rigorous stability analysis is beyond the scope of this paper, we will argue that 
a finite energy barrier will always separate the two surface-bound states. This involves the construction 
of a natural homotopy connecting these states. Let 

Y t = (l-t)Y + (l-*)R(o;)Y lj (59) 

where we represent by Y the initial state with n = 2,p = 1 and by Yi the final state with n = 3,p = 2. 
These interpolating non-equilibrium states will also be confined within the sphere due to the convexity 
of the latter. The two boundary states have R = 2{^] A rotation R(uj) of the final state about the 
axis of symmetry is introduced in order to initialize the pointwise sum so as to avoid the development 
of cusps. Intermediates in this homotopic sequence are illustrated in Fig. [8] with the constant choice 
uj = 7r/6. The bending energy is plotted in Fig. [9] as a function oft. It exhibits a finite potential barrier 
separating the two states. While this does not prove that the n = 2 state is stable in this regime, it 
does suggest that it is. 

One cannot rule out the existence of semiattached equilibrium states of the loop. We were, however, 
unable to construct any such states. If they do exist, one would expect them to be unstable. 



5.4 Transmitted forces 



Using Eq. (51), the local force per unit length transmitted to the sphere is given explicitly as the 
following function of s, 

(60) 



X(s) = ^ 2 g + a = q 2 (2dn 2 (qs,m) - l) , 



where dn 2 [u,m] = 1 — msn 2 [u,m] [24] . Whereas the transmitted force depends on the local value of 
k 9 , both Kg and a depend on the boundary conditions associated with closure. Thus A does depend 
on the global loop geometry. Note that the expansion of the expression in (60) to second order in 
AR = R — 1 reproduces the result for weak confinement (Eq. (|46|)). 



In Fig. 10 the maximum and minimum values of A have been plotted as a function of R for the two- 
fold "ground" state and its descendants. It is positive everywhere with the minimum given by a. One 



17 As written down, the homotopy does not preserve length. This can be achieved by stretching Y and Yi appropriately 
at intermediate values of t. 
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1.0 1.5 2.0 2.5 3.0 3.5 4.0 

R 



Figure 7: Comparing the geodesic energies of the states with two and three fold symmetries in an interval of 
R in the neighborhood of R = 2 (solid and dashed lines respectively). The geodesic energy, and thus the total 
energy, of the three-fold is lower than that of the two-fold in a finite interval of R beginning at R = 2. 




(d) t = 0.6 



(e) t = 0.8 



(f)t=l 



Figure 8: (Color online) Homotopy connecting the two states n = 2,p = 1 and n = 3,p = 2 with R = 2. 
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Figure 9: Total energy along a homotopy interpolating between two stable ground states. An energy barrier 
always exists between the two. 

observes that A does not behave in a monotonic way, suffering a positive discontinuity whenever R = p, 
where p is an odd integer. This is associated with the transition from orbital to oscillatory behavior at 
these values and is the analog of the Euler instability described in the discussion of a weakly confined 
loop. If the unstretchability constraint is relaxed the discontinuity will be smoothed but the jump will 
persist. Within the individual intervals between discontinuities both the maximum and minimum values 
decrease monotonically with R; this is associated with the reduction in the force once buckling into a 
state of oscillation has occurred. Curiously, whereas the maxima within these intervals decrease with R, 
the corresponding minima increase. This occurs in such a way that the mean of their values increases 
monotonically. Both maxima and minima approach the value a = 1 asymptotically. 



Recall that the transmitted force is bounded from below by a (see Eq.(45)). Thus, if a > 0, then A is 
also; the confined loop will then push on the sphere everywhere. If, however, a < then A may change 
sign along the curve. While a is positive in the ground state and its descendants, it may become negative 
in the excited states of sufficiently long loops. In particular, we find that all states with n > 5 and p = 1 
exhibit regions in which A turns negative at values of R below the onset of self-contact. When n = 5 
and p = 1, this occurs when R > 2.536 (see fig[TT](a)). Such states are, however, unstable. It may 
appear counterintuitive that A can be negative. It has to be remembered, however, that while A depends 
only on the local geometry, this geometry is determined by the global behavior of the loop through the 
boundary conditions. If one cuts a loop somewhere, the state will immediately relax into a geodesic 
state with constant positive A everywhere. The existence of regions along the loop where A < does 
not necessarily signal a tendency to detach. This will depend on details of the energy landscape. 



The states examined here are those consistent with the bound Eq. (23). As was pointed out, states 
that violate this bound are always unstable with respect to the unbinding of the loop into the interior. 
Such a state is illustrated in Fig. pT|(b), making two self-intersecting orbits in the northern hemisphere. 
The curvature is high within the two orbits so that this state is evidently an excited state of the loop. 
The transmitted force A is negative everywhere below the circle with k 9 = \/—2g (see Fig. [TT|(b)). 
This suggests that the likely mode of instability will involve the loop unbinding along regions where A 
is negative, allowing the loop to relax toward the ground state n — 2, p = 1 through the unwinding of 
the high curvature orbits. 
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Figure 10: Transmitted force. Maximum (solid line) and minimum (dashed line) values of A as a function of 
AR in the two-fold ground state with n — 2, p = 1, A in the two-fold ground state and its minimum energy 
descendants. 

Intuitively, one might have expected the total force transmitted to the sphere, A = J dsX to increase 
with loop size. This is not generally the case, except in the limit of large loops, where A increases linearly 
with R, a consequence of the asymptotic behavior of a. Notice also that A is not the same as the energy 
H (divided by the radius of the sphere) for finite values of R. One has A = H — 2nR(l — c(R)). In 
the limit, however, the two do coincide. A is plotted as a function of R in Fig. [12} 

6 Discussion 

The confinement of a closed semiflexible loop by a surface has been examined. The constraint on 
the spatial degrees of freedom of the loop is enforced in the variational principle using the method of 
Lagrange multipliers. The loss of Euclidean invariance of the energy of the unconstrained loop under 
confinement is quantified by the Lagrange multiplier which, in equilibrium, gets identified as the force 
transmitted to the surface by the confined loop. 

We have focused on a description of the ground state of a confined loop. If the loop is short this consists 
of an oscillation about a geodesic circle, exhibiting a two-fold dihedral symmetry. This is the only stable 
state of the loop in this regime. The description of the ground state and its excitations gets more 
complicated when the loop length is increased. When R is increased above 2R 0l a new set of states 
comes into existence that oscillates about a double covering of a geodesic circle. Among these the state 
with three-fold symmetry has lowest energy; within a finite band of values of R } this state also replaces 
the two-fold as the ground state. The two states remain separated by an energy barrier so that the 
two-fold remains stable. As the loop size increases these states get replaced by descendants with higher 
n-fold symmetries which alternate as ground states. Both the energy and the transmitted force suffer 
discontinuities associated with changes of symmetry. They do not generally increase monotonically with 
loop size except asymptotically. In this limit the energy gap between the ground state and excited states 
disappears. 

Our examination of the confinement of a loop within a spherical cavity might lead one to expect that 
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Figure 11: (Color online) (a) The state n = 5, p = 1 with R = 2.536. Above this value of R, A may be 
negative in places, (b) The state n = 2, p = 1 with R = 1.85 inhabits a single hemisphere. These orbital states 
always have intervals where A is negative. 

the equilibrium states of a confined loop will always attach. However, other confining geometries display 
very different behavior. Consider, for example, the confinement of a loop by a cylinder of smaller radius. 
Analogs of the n-folds exist but, in general, they do not provide the ground state as one can easily 
verify by playing with metal rings in a wastepaper basket. The details are surprisingly complicated. In 
this context, it would be interesting to understand how the spherical equilibrium states "evolve" under 
elliptical deformations of the sphere, or what effect surface irregularities or pores have on the confinement 
process. 

When contact with the confining surface is incomplete, one needs to address the issue of boundary 
conditions at points of contact. If the contact is due entirely to geometric hindrance, the boundary 
conditions at isolated points of contact are simple: the tangent vector to the curve is tangent to the 
surface at these points. If contact extends over a finite region, however, in addition to the tangent 
vectors the curvatures will also be continuous at the boundary of the region of contact. There may 
also be a tendency to either promote or inhibit adhesion to the confining surface. A simple way to 
accommodate such interactions is to introduce a local contact energy, proportional to the length of 
the contact interval. The boundary conditions will reflect this additional energy. While the curvature 
will suffer a discontinuity analogous to the discontinuity associated with the contact line bounding the 
region of contact between a fluid membrane and an attractive substrate [26,12711281129], this is not the 
full story. The normals may rotate about the tangent vector, aligning themselves with some preferred 
direction on the confining geometry. The extension of our framework to accommodate contact energies 
will be presented elsewhere. 

In this paper the role played by the confining surface has been passive: we have not considered the 
possibility that it may deform in reaction to the presence of the confined polymer. Surface deforma- 
tions due to membrane bound polymers may also play a role in shaping the morphologies of biological 
membranes. The forces transmitted to the membrane will now provide a source for the surface stress. 
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Figure 12: A vs. R (solid line). The dashed line indicates the contribution from geodesic curvature A g . The 
dotted line is the contribution = 2nRa{R). 

Understanding this coupling poses technical challenges. What is clear is that there will be interesting 
physics associated with the competing elastic energies [30] . 
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Appendix A Tension in an unconstrained elastic space curve 



There are various derivations of Eq. (19). The approach adopted here involves treating the tangent 
vector to the curve T = Y' as an independent variable. If the curve is also parametrized by arc-length, 
T is a unit vector, T 2 = 1. Thus, consider the energy functional (k 2 = T /2 ), 

ff 1 [Y,T,A,F] = J da Q T' 2 - ^A(T 2 - 1) + F • (T - Y')) ; (A.l) 

the presence of the two constraints liberates T to be varied independently. Now the variation of Hi 
with respect toY is given by 

SyHx = ! ds F' • 5Y , (A.2) 



(A.3) 



whereas its counterpart with respect to T is 

St Hi = J ds (F - T" - AT) • 5T . 



Modulo the Frenet-Serret equations for the curve, S T Hi = in Eq. (A.3) implies 

F = (-K 2 + A)T + k'N + ktB. 



(A.4) 
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For an isolated elastic curve, F' = 0. If the curve is constrained to lie on a surface, it was seen in Sec. 
[2] that F' = An. Thus the tangential equation t • F' = continues to hold. This equation can be 
integrated to determine the Lagrange multiplier A associated with the parametrization by arclength 



A 2 

A -k z -c. 
2 

where c is a constant of integration. Equation ([4]) follows. 
With respect to the Darboux frame, F is given by 



(A.5) 



F = 



f 9 2 — C\T + (< + KgTg) II + (K' g - K n Tg) 1 . 



(A.6) 



where we have used the relations (B.4) derived below. 



Appendix B Darboux frame for surface curves 

The structure equations (analogous to the Frenet-Serret equations) describing the Darboux frame are 
given by 

T' = K n n + K g l, (B.la) 

n' = -K n T-r g l, (B.lb) 

I 7 = -K g T + TgIL. (B.lc) 

The normal curvature, the geodesic curvature, and torsion are given, respectively, by 

k u = r • n = -K ab t a t b , K g = r • 1 = rt b v h t a , Tg = r • n = -i^ a / 6 . (b.2) 

Here t a and / a are the components of the vectors T and 1 with respect to a basis of surface tangent 
vectors adapted to the parametrization, e ai a = 1, 2: T = t a e ai 1 = l a e a . K ab = e a • c^n is the extrinsic 
curvature tensor defined on M and V a is the covariant derivative compatible with the induced surface 
metric g ab = e a • e b . Whereas K n and r g depend on the extrinsic curvature, k q is defined intrinsically; it 
depends only on the surface metric g ab . The Frenet curvatures are related to their Darboux counterparts 
by 

Kg — k sin a; , K n = k coscj , r = uJ — r g , (B-3) 

so that k 2 = + Thus, the Frenet curvature can be decomposed into its intrinsic and extrinsic 
parts. The Frenet torsion r is the sum of the derivative of the angle connecting both frames and the 
geodesic torsion. Note that r g involves two derivatives, whereas r involves three. The extra derivative 
is associated with uj f , the rotation rate of one frame with respect to the other. 



The identities Eqs. (B.3) imply the following: 



Kg = —Kg + K n (T + Tg) , K^ = —K n - Kg(r + Tg) . (B.4) 

K K 

These relations are used in Sec. [2] to express the Euler-Lagrange derivatives for the curve in terms of 
Kg, K ni and r g and their derivatives. 
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Appendix C Identities for weak spherical confinement 

We first derive the relationship between arc-length and the angle (p correct to second order. One has 



s — I dip 



dip J 



+ sirr i? 



1/2 



Using Eq. (29) and the harmonic approximation for k 9 given by Eq. (38), one obtains 



sin d 



MS 



cos rap 



dt) 1 dn\ nA\ 
and — ~ — — - — = — — — sin no? . 

dip M dip M 



Thus, in this approximation, 



s ~ <p + 



1 A\ 



4 [n 2 - l) 2 



(n 2 - l)p 



(n 2 + l) . 



2n 



sin2n^ 



This implies thatf^ 



AR 



1 



Ai 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



4(n 2 - 1) 

We now derive Eqs. (40) and (41). To do this one needs to examine the boundary conditions p4| ) and 
(35) correct to second order in e = V AR. At this order, the turning point k\ of the potential V(k 9 ), 
defined by Eq. (25), coincides with A\ given in the harmonic approximation by Eq.(39). Thus, the roots 
in the factorization of E 2 — V{k) given in Eq. (26), where E is is the constant defined below Eq. (25) 
are 



k\ = 2{E -<t)bsA\, K 2 = 2(E + a) 



-An 1 



(a 2 - M 2 ) - 2(a 2 + M 2 ) . 



One thus has, correct to second order, 
1 1 1 



V E2 - V(K g ) 



n 



1 - ^r(^ 2 - M 2 ) - -^(a 2 + M 2 ) 



4n 4 



4n 2 



8n 2 



Equation (34) then implies that 



Al 



Furthermore, to second order, 

\K 2 g + (J-l 



M 



M 2 - k\ 



1 + 



n 2 — 1 



fcr 2 - M 2 ) + 



n 2 + 1 



2(n 2 -l)n 2 -l 9 



Equation (35) then implies that for 1, 

1 



((7 2 - M 2 ) 



1 n 2 + 1 2 



4(n 2 - 1) n 2 - 1 



(C.5) 



(C.6) 



(C.7) 



(C.8) 



(C.9) 



18 Ai vanishes when n = 1, a solution which is identified as a trivial rotation of the equatorial loop about an axis in 
the equatorial plane. 
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